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Abstract. In recent years, the Lyman-a absorption observed in the spectra of high-redshift 
quasars has been used as a tracer of large-scale structure by means of the three-dimensional 
Lyman-a forest auto-correlation function at redshift z ~ 2.3, but the need to fit the quasar 
continuum in every absorption spectrum introduces a broadband distortion that is difficult to 
correct and causes a systematic error for measuring any broadband properties. We describe a 
k-space model for this broadband distortion based on a multiplicative correction to the power 
spectrum of the transmitted flux fraction that suppresses power on scales corresponding to 
the typical length of a Lyman-a forest spectrum. Implementing the distortion model in fits for 
the baryon acoustic oscillation (BAO) peak position in the Lyman-a forest auto-correlation, 
we find that the fitting method recovers the input values of the linear bias parameter hp and 
the redshift-space distortion parameter Pf for mock data sets with a systematic error of less 
than 0.5%. Applied to the auto-correlation measured for BOSS Data Release 11, our method 


improves on the previous treatment of broadband distortions in BAO fitting by providing a 
better fit to the data using fewer parameters and reducing the statistical errors on j3F and 
the combination 6_p(l + j3F) by more than a factor of seven. The measured values at redshift 

^ — 2.0 are pp — J-•oy_QJQ _o.i9 -0.28 + PF ) — ’-'•'^*^-0.007 -0.014 -0.022 

and 3a statistical errors). Our htting software and the input files needed to reproduce our 

main results are publicly available. 

Keywords: dark energy experiments, Lyman alpha forest, cosmological parameters from 
LSS, baryon acoustic oscillations 
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1 Introduction 

Baryon acoustic oscillations (BAO) is one of the most well-established observational probes 
of cosmic acceleration [1]. The characteristic length scale of BAO imprinted in the large- 
scale clustering of the density field provides a standard ruler that can be used to study the 
expansion rate and geometry of the Universe as a function of redshift [2]. A measurement of 
the BAO scale in the radial direction probes the Hubble parameter, whereas a measurement 
in the transverse direction provides constraints on the angular diameter distance. Until 
recently, the BAO technique had mainly been exploited using the clustering of galaxies at 
redshift z < 1 [3-7]. The Baryon Acoustic Oscillation Survey [8] (BOSS) of the Sloan Digital 
Sky Survey III [9, 10] (SDSS-III) reported a 1% measurement of the cosmic distance scale 
at 2 : = 0.57 and a 2% measurement at z = 0.32 [11] using a sample of nearly one million 
galaxies taken from SDSS-III Data Release 11 [12] (DRll). 

The expansive BOSS data set includes an unprecedented number of high-redshift quasar 
spectra that enable us to explore the large-scale structure at z > 2 through the Lyman-a 
(Lya) forest technique. The spatially correlated absorption of Lya photons by the neutral 
hydrogen component of the intergalactic medium along the lines of sight to quasars traces the 
matter distribution, thus allowing for measurements of BAO during the epoch of decelerated 
expansion when the universe was strongly matter dominated [13, 14]. The first observation 
of three-dimensional large-scale correlations in the fluctuations of the Lya forest absorption 
using BOSS quasar spectra was reported in [15], providing observational constraints on the 
Lya forest linear bias parameter and redshift-space distortion parameter at z = 2.25. Using 
a sample of approximately 50,000 BOSS quasar spectra, the first detection of the BAO peak 
feature in the Lya forest auto-correlation function at z ~ 2.3 was presented in [16-18]. The 
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analysis was improved to yield a more precise measurement of the BAO peak position using 
the three-times larger DRll quasar sample [19]. Complementary measurements of BAO at a 
similar redshift can be achieved by considering the correlated clustering of the intergalactic 
medium with the quasars themselves. The three-dimensional large-scale cross-correlation of 
quasars and the Lya absorption was detected in [20] and subsequently extended to larger 
separations for DRll to measure the BAO scale [21]. By combining the BAO results for the 
DRll auto- and cross-correlation, the Hubble parameter and the angular diameter distance at 
z = 2.34 were measured to 2% and 3% accuracy, respectively. The cosmological implications 
of the BAO measurements from the BOSS galaxy clustering and Lya forest analyses, in 
combination with data from other probes, were studied in [22]. 

A method for fitting Lya forest BAO was developed in [18]. Here, we refine the method 
by constructing a model for the distortions of the broadband shape of the measured cor¬ 
relation function introduced by the analysis. Broadband distortions primarily arise due to 
continuum fitting in the Lya forest of quasar spectra. The standard approach for accom¬ 
modating these distortions in BAO fitting is to introduce a smooth broadband function in 
r-space that is additive and/or multiplicative on the correlation function. The broadband 
function is typically required to feature many free parameters (usually about ten) in order 
to provide a good fit to the measured correlation function. In the method developed here, 
the broadband distortion is instead modeled in A:-space as a multiplicative correction to the 
power spectrum of the transmitted flux fraction that suppresses power on scales correspond¬ 
ing to the typical length of the Lya forest. This distortion model has a clearer physical 
interpretation and allows for a much more precise determination of the Lya forest linear bias 
parameter and redshift-space distortion parameter. This approach additionally provides a 
better fit to the Lya auto-correlation using fewer parameters. 

Our analysis uses the DRll Lya auto-correlation from [19] based on the BOSS spectra 
of 137,562 quasars in the redshift range 2.1 < z < 3.5. The BAO fit results reported in 
that paper can be directly compared to the results obtained here using the distortion model, 
because changes are only made in the BAO modeling and not in the data. The fit input files 
and the fitting code are made publicly available, so that readers may reproduce our main 
results. 

This paper is organized as follows: Section 2 describes the BAO model, including models 
of the non-linear correction of the transmission power spectrum. In section 3, we define our 
model of the broadband distortion due to continuum fitting. Section 4 explains the method 
for fitting the BAO model to the measured auto-correlation function. In section 5, we test the 
distortion model using auto-correlation functions for mock data sets, and section 6 describes 
the fits to the DRll auto-correlation and presents our main results. The paper concludes 
with a discussion in section 7. 

2 BAO model 

In this section, we describe our method to calculate a theoretical prediction of the Lya auto¬ 
correlation that can be fitted to the measured auto-correlation to determine the BAO scale. 
The modeling is done in fe-space by means of the Lya forest transmission power spectrum, 
which is Fourier transformed into the Lya auto-correlation. Our fiducial cosmology is the 
same as in [19], namely a spatially flat ACDM universe with Dm = 0.27, h = 0.7, = 

0.0227, = 0.00064 (massive neutrinos), erg = 0.79 and Ug = 0.97. 


- 2 - 


2.1 Transmission power spectrnm 

We model the power spectrum of the Lya forest transmitted flux fraction at redshift z, 
including linear redshift-space distortions, non-linear effects and broadband distortion due 
to continuum fitting, as 

PF{k,Hk,z) = b%{z) [l + PNhik, Hk, fj,k)Dcik\\) , (2.1) 

where k is the modulus of the wavenumber, is the cosine of the angle of the wavenumber 
from the line of sight, and k\\ = k^k is the component along the line of sight. The model builds 
on the linear theory prediction (e.g., [13]), for which Pf is reduced to the redshift-dependent 
bias factor bF{z) of the Lya transmission, the term depending on the redshift-space distortion 
factor 13F (which in principle depends also on redshift, but is assumed to be constant in this 
paper), and the linear matter power spectrum. Here, we separate the non-linear effects 
modifying the power spectrum into a multiplicative non-linear correction for small scales, 
D-^h{k, that is described in section 2.2, and a smoothing term that models the effect of 
large-scale bulk motions on the BAO peak described below. The term Ilc(^||) rnodels the 
broadband distortion due to continuum htting, and is described in detail in section 3. 

The non-linear matter power spectrum is modeled by 

Pi^hik, Hk, z) = exp [-A;^S^(/Xfc)/2] Pi,{k, z) , (2.2) 

where Ph{k, z) is the (isotropic) linear matter power spectrum obtained from CAMB [23] for 
our fiducial cosmology. Non-linear effects on the matter density field due to large-scale bulk 
velocity flows are modeled with an anisotropic Gaussian roll-off of the linear matter power 
spectrum [24], with 

S^(/rfc) = PkP'l + (1 “ ^1)^1 ■ (2-3) 

The effect of the anisotropic non-linear broadening is a convolution of the correlation function 
that broadens the BAO peak. The parameters quantifying the strength of the radial and 
transverse broadening are related by 


^ = 1 + /, (2.4) 

where / = d{ln g) /d(ln a) « 0[(;^^(z) [25] and g{z) is the growth factor. We adopt the 
values Sji = 6.41 Mpc//i and S_l = 3.26 Mpc/h expected at z = 2.4 [18]. The non-linear 
broadening model provides a valid description of the effect on scales of the BAO feature, 
but applies unphysical filtering of small-scale structure. A more refined model would feature 
scale dependent broadening parameters. 

The redshift evolution of the linear matter power spectrum at z > 2 (where / ~ 1) is 
governed by the growth factor with respect to a reference redshift Zj-ef such that 

P^ik, z) = P^ik, Zrei)g‘^{z) = Phik, Zref) ' ( 2 '^) 

Following the result of [26], we assume that Pf evolves with redshift according to the com¬ 
bination ^ ^ 

h\{z)g^{z) = b\{zrei) ) • ( 2 . 6 ) 

+ Zrei ) 
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Combining equation (2.5) and equation (2.6) produces the result that the linear bias evolves 
as ^ g 

6i.(z)=6i.(zref)f:r^^V . (2.7) 

\ i + Zref J 

We will use the shorthand notation hp = bpizrei) in the remainder of this paper. 

The three-dimensional Lya forest auto-correlation function ^p{f,z) is most generally 
related to Pp{k, z) by a Fourier transform 

^F{f,z) = j d^kPp{k, z) exp{-ik ■ r) , (2.8) 

which can be evaluated numerically with a three-dimensional Fast Fourier Transform (3D 
FFT). Specializing to power spectra that are cylindrically symmetric, and therefore only 
functions of A:|| and k±, the polar angle integration can be performed explicitly, 

1 

^Fir\\,r±,z) = — dk±k±Jo{k±rj_)F{r\\,k_i,z) , (2.9) 

27r Jq 

where F{r^^,k±, z) is the ID Fourier transform of Pp along k^^ only 

2 r+oo 

F{r\\,k±,z) = — J dk\\Pp{k\\,k±,z)exp{-ik\\r\\) . (2.10) 

This result suggests an alternative two-step approach to the numerical computation; first, 
perform N^^N± ID FFTs, followed by A/^||A^_l separate ID numerical integrations, where 
and A^_l are the number of grid points along the radial and transverse axis, respectively. 
We are currently using the 3D FFT method, for which we extract ^p{r^^,r±, z), and leave 
a two-step implementation for future work. Our 3D FFT uses grid spacing Ar = 4 Mpc/h 
and N = 400 grid points along each axis, corresponding to coverage up to fcmax = TrjAr = 
0.79 /i/Mpc with spacing Ak = 2 p/{NA r) = 3.9 x 10“^ /i/Mpc. This resolution is sufficient 
to reach a satisfactory accuracy of the correlation function calculation when fitting for the 
BAO peak position. 


2.2 Non-linear correction 


The function Dnl in equation (2.1) represents a correction due to non-linear effects that 
are dependent on the non-linear clustering of matter, and the hydrodynamic and thermal 
evolution of intergalactic gas. These non-linear effects can only be properly modeled through 
hydrodynamic simulations of the intergalactic medium. A parameterized formula for the non¬ 
linear correction DNL(fe) fJ-k) can then be constructed by fitting to the ratio of the calculated 
Pp for the simulation to Pp. We note here that the anisotropic non-linear broadening that 
we include in the Gaussian term in equation (2.2) arises from fluctuations at the BAO scale, 
which are not captured in the simulations that have been used so far for calculating Dnl 
because of their small box size. We can therefore safely assume that the non-linear broadening 
is a separate effect from anything that is already included in Dpip, although this will not be 
true with future simulations with larger box size. In this paper, we employ two non-linear 
correction models from the literature for separate parts of the analysis. 

The first model, which we designate DNLl, is the fitting formula of [13], 


Dj^p{k,Hk) = exp 



Pnl 



kfj-k 

ky{k)) 


( 2 . 11 ) 


- 4 - 







z 

5NL 

ky 

ay 

by 

kp 

3.0 

0.792 

1.16 

0.578 

1.63 

17.1 

2.8 

0.773 

1.16 

0.608 

1.65 

19.1 

2.6 

0.781 

1.15 

0.611 

1.64 

21.0 

2.4 

0.851 

1.06 

0.548 

1.61 

19.5 

2.2 

0.867 

1.05 

0.514 

1.60 

19.3 


Table 1. Parameter values as a function of redshift for the non-linear correction model DNL2. The 
parameters fcy and fcp are given in units of /i/Mpc. 


where 

/ /c 

kyik) = kyo ^1 -h p-J . (2-12) 

The DNLl model has eight parameters and we use the fixed values from table 1 of [13] 
derived for b'p = 0.0173 and Pf = 1-58 at z = 2.25: k^L = 6.40 /i/Mpc, pnl = 0.569, 
kp = 15.3 h/Mpc, pp = 2.01, kyo = 1.220 /i/Mpc, py = 1.50, ky = 0.923 /i/Mpc and 
Py = 0.451. The DNLl model is used exclusively in our tests on mock data sets, as explained 
in section 5. Since these parameter values are consistently used for both generating and 
analyzing the mock data sets, it makes no difference that they were derived for a different 
fiducial cosmology. 

The second model, denoted DNL2, is the fitting formula of [27], 


DML(k,iUk) = exp 


gNL^^(k) I 1 - 


_k_ 

ky 


av 




_k_ 

kp 


where 


A\k) = 


k^PUk) 

27r^ 


(2.13) 


(2.14) 


The DNL2 model has five parameters and depends on the linear matter power spectrum 
Pp{k). We use the DNL2 model in our fits to the DRll auto-correlation in section 6. Table 1 
lists the parameters ^nl, ky, ay, by and kp as a function of redshift; we interpolate linearly 
in each parameter to the effective redshift of the BAO measurement. The simulations of 

[27] feature values of the cosmological parameters that are consistent with Planck results 

[28] , which makes transferring the derived DNL2 parameters to our fiducial cosmology non- 
optimal. However, the differences are small once the amplitude of the linear matter power 
spectrum is re-scaled by the factor (0.834/0.79)^ to match their value as = 0.834. 

Both model parameterizations feature three terms inside the exponential function which 
represent different physical effects: the isotropic enhancement of power due to non-linear 
growth (subscript “NL”); the isotropic suppression of power due to gas pressure below the 
Jeans scale (subscript “P”); and the suppression of power due to line-of-sight non-linear 
peculiar velocity and thermal broadening (subscript “V”). 

Figure 1 illustrates the non-linear correction models implemented for analyzing the 
mock data sets (DNLl) and the DRll data (DNL2). DNLl is fixed at z = 2.25, while DNL2 
is shown for z = 2.3, closer to the effective redshift of the DRll BAO measurement. The 
fact that there are some differences between the models is not surprising, because they were 
derived using somewhat different techniques for simulations with different properties and 
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k (h/Mpc) 


Figure 1. Illustration of the non-linear correction models DNLl (green) and DNL2 (orange). The 
solid, dashed and dotted lines show the models for = 0.9, fj.^ = 0.5 and = 0.1, respectively. 
DNLl adopts the parameter values at z = 2.25 specified in the text. For DNL2, we have used 
parameter values at z = 2.3: (/nl = 0.859, ky = 1.055 /i/Mpc, ay = 0.531, by = 1.605, fcp = 
19.4 h/Mpc. 


fiducial cosmologies. The most important qualitative difference between the models is that 
DNL2 quickly converges to unity for decreasing k, as one would expect in the linear regime, 
whereas DNLl converges more slowly and thus produces a non-linear contribution even on 
large scales. In our application, the 3D FFT resolution used for the BAO fits restricts the 
wavenumber range to A; < 0.79 /i/Mpc. 

Figure 2 shows the effects of the non-linear correction and the anisotropic non-linear 
broadening on Pp and the corresponding correlation function at z = 2.3. The increase 
of power on large and intermediate scales by DNLl affects the broadband shape of the 
correlation function. In comparison, DNL2 produces a modest change of the shape. The 
suppression of power on small scales by the non-linear broadening smooths the correlation 
function to modify the broadband shape and broaden the peak. 

3 Distortion model 

Fitting the continuum level in the Lya forest of quasar spectra is an operation that induces 
distortions of the broadband shape of the correlation function. Our analysis is based on 
the DRll auto-correlation for which the Lya forest continua were fitted using the “C2” 
method [19]. In this method, the Lya forest continuum Cq for each quasar is modeled as 
the product of the mean continuum C for the quasar sample, determined by stacking the 
normalized spectra in the rest frame, and a wavelength-dependent function that sets the 
overall continuum amplitude and slope, 

Cq{X) = [aq + bqlog{X)]Ci\r{) . (3.1) 

The parameters aq and bq are fitted using a maximum likelihood approach that imposes that 
the transmitted flux fraction F, calculated as the ratio of the observed flux with the con¬ 
tinuum, follows an assumed log-normal probability distribution function [16]. The spatially 
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Figure 2. Effects on the fc-weighted power spectrum Pp (left-hand panels) and corresponding r^- 
weighted correlation function (right-hand panels) from the non-linear contributions at z = 2.3. 
Top panels show the results for /i = 0.9 and bottom panels for /r = 0.5. Curves indicate the results 
for no non-linear effects (solid, black), defined by Il|| = S_l = 0 and Dnl = 1; only including the 
anisotropic non-linear broadening (dotted, black); only including the non-linear correction model 
DNLl (solid, green); only including the non-linear correction model DNL2 (dashed, orange). We have 
adopted j3p = 1.4 and bp{l + j3p) = —0.351 (corresponding to hp{l + l3p) = —0.336 at z = 2.25 [15]). 
Distortion due to continuum fitting has not been included [Dq = 1). 


correlated fluctuations in the transmission field are calculated as 

where the mean transmitted flux fraction F is determined by calculating the average of all F 
values at each wavelength, such that, by construction, the mean value of hi? vanishes at each 
wavelength. Broadband distortions arise because the continuum fitting cannot distinguish 
between a change in the continuum amplitude and a true offset in and will typically favor 
the former to adjust out the large-scale Fourier modes of <5ir. The exact characteristics of the 
broadband distortion depend on the details of the continuum fitting method as well as the 
properties of the quasar sample, meaning that the distortion will be different if one changes 
the fitting method or modifies the sample. 

Analytically modeling the broadband distortion introduced by continuum fitting is com¬ 
plicated, because a large number of quasars with different properties and over a wide redshift 
range are used to measure the correlation function. There can be intrinsic variations in the 
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quasar spectral shape and flux calibration systematic errors, and the distortion depends on 
how these variations are allowed in the method to fit the continuum. In addition, the phys¬ 
ical length of the Lya forest changes with redshift, with weights at each spectral pixel that 
vary depending on each particular quasar observation. The approach that has been taken so 
far is to use analytical prescriptions for estimating an additive distortion correction to the 
correlation function [15, 29], but this has suffered from modeling uncertainties that degrade 
the accuracy to measure large-scale, broadband properties of the correlation function. 

In this paper, we adopt a different approach based on noticing two facts about the 
distortion correction. First, in /c-space, the correction should be a multiplicative function 
that depends only on the radial component fey. Second, this function should suppress the 
large-scale power on the scale that corresponds to the typical length of the Lya forest in 
our quasar spectra. To derive the first property of the distortion correction, we consider the 
impact of this distortion on single Fourier modes 5F{k\\-,k±) of a fixed k\\ but different k_\_, 
approximating the quasar spectra to be along nearly parallel lines of sight. The variations 
in A;_i_ are only affecting each individual spectrum by a phase shift, but this phase shift is 
uniformly distributed from 0 to 27r independently of the value of k_\_. Since the continuum 
is determined for each individual spectrum without any reference to nearby parallel spectra, 
the average suppression of the power can only depend on k\\. The correction is multiplicative 
since a change in the undistorted model (e.g., a change in bp) should scale correspondingly 
in the distorted model. The second property is a consequence of fitting the continuum on 
each spectrum over the length in which the Lya forest absorption is observed, which sets the 
characteristic scale determining which Fourier modes are removed most effectively. 

The DRll auto-correlation is based on the spectra of 137,562 quasars in the redshift 
range 2.1 < z < 3.5. Using the rest-frame wavelength range 1040 < Arf < 1200 A as the width 
of the Lya forest, the maximum forest length (i.e. the comoving separation spanned by the 
first and last spectral pixel of the forest) is approximately 635 Mpc//i at z = 2.47. The length 
decreases at lower redshift, because the spectra are truncated at observed wavelength A = 
3600 A near the edge of the BOSS spectrograph [30]. At higher redshift, the length decreases 
slowly because of the increasing Hubble parameter. Figure 3 displays the distribution of the 
forest lengths for the quasar sample, with a mean value of approximately 520 Mpc/h. The 
scale of the typical forest length corresponds in A;-space to k\\ ~ 0.01 /i/Mpc, so we expect 
our distortion function to vary on this characteristic scale. 

We construct the distortion model Dc{k\\) as a sigmoid-type function that is zero at 
/cji = 0 and smoothly increases to approach unity on small scales. We consider two parame- 
terizations and compare their performances using fits to auto-correlation functions for mock 
data sets in section 5.3. The first model, denoted DCl, is given by 


Dcik\\) = tanh 



and the second model, designated DC2, is 


Dc{k\\) 


(fc||/fcc + l)"^' 


1 


PC 


{k\\/kc + l) + 1 


(3.3) 


(3.4) 


Both distortion models have two free parameters, fee and pc, that are fitted simultaneously 
with the parameters of the BAO model. Because of the differences in parameterizations, 
the meanings of fee and pc are not exactly the same between the models. As described in 
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Figure 3. Distribution of the Lya forest lengths for the BOSS DRll quasars in the redshift range 
2.1 < 2 < 3.5 that were used to measure the auto-correlation. The adopted bin size is 20 Mpc/h. The 
calculated forest lengths are the maximum comoving separations given the quasar redshift and the 
forest rest-frame wavelength range 1040 < Arf < 1200 A. The wavelength cut at observed wavelength 
A = 3600 A, near the edge of the spectrograph, truncates the spectra of quasars with z < 2.48. 


section 5.3, we fix the pc parameter for DC2 based on the fit results for mock data sets, in 
order to neutralize its degeneracy with kc- 

Figure 4 illustrates the distortion models for varying parameter values. DCl shows a 
clear distinction in the effect of the parameters, where kc controls the scale of the suppression 
and Pc modifies the shape. In the case of DC2, the difference between varying kc and pc is 
more subtle on the scale of the suppression, indicating a degeneracy between the parameters. 
Figure 5 shows the effects of the distortion model DC2 on Pp and the corresponding correla¬ 
tion function at z = 2.3. The suppression of the power spectrum on large scales distorts the 
broadband shape of the correlation function, but does not change the position of the BAO 
peak. 

4 BAO fitting method 

We implement the models for the non-linear correction and the broadband distortion into the 
publicly available fitting software baofit^. Our method for fitting the BAO in the Lya auto¬ 
correlation is largely the same as in previous analyses using BOSS data [16-19, 21]. The 
correlation function is decoupled into a smooth and a peak component, 

Cir(r||,r_L,a||,a_L) = 4mooth(r||, r_L) -b Opeak • ^peak(a||r||, a_Lr_L) . (4.1) 

Decoupling the peak is necessary in BAO fitting when the broadband shape is insufficiently 
understood. Although the fc-space distortion model developed in this paper greatly improves 
the modeling of broadband distortions from continuum fitting for a decoupled fit, it is not 
quite refined enough for an unbiased measurement of the BAO scale for a fit where the peak 
is coupled to the broadband component. The parameter Upeak controls the amplitude of the 

^See Appendix A for information on where to download the fitting software. 
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Figure 4. Illustration of the distortion models. Left: Distortion model DCl for parameter values 
kc = 0.02 h/Mpc and pc = 1 (solid, red). The gray and red bands are the regions enclosed by 
changing the kc and pc parameters by ±50% (the gray and red dashed lines indicate the larger 
values). Right: Distortion model DC2 for parameter values kc = 0.003 h/Mpc and pc = 3 (solid, 
blue). The gray and blue bands are the regions enclosed by changing the kc and pc parameters by 
±50% (the gray and blue dashed lines indicate the larger values). 


BAO peak and is kept fixed to its nominal value Opeak = 1- The anisotropic shift of the 
observed peak position relative to the fiducial peak position is described by the radial and 
transverse scale factors 

^ [DHizes)/rd] ^ [DAizes)/rd] , , 

where Zes is the effective redshift of the BAO measurement, Dh{z) = c/H{z) is the Hubble 
distance, Da{z) is the comoving angular diameter distance, and is the sound horizon at 
the drag epoch [31] with a measured value 147.49 Mpc reported by Planck [28]. 

Decoupling the correlation function into a smooth and a peak component for the BAO 
fitting means that we must perform the corresponding split in the model of the matter power 
spectrum. Because the anisotropic non-linear broadening model is mainly valid on BAO 
scales, we choose to apply the broadening to the peak only, such that 

-Pnl(^, Aifc, z) = exp(-/c2S^(/rfc)/2)Ppeak(fc, z) ± Psmooth(A:, z) . (4.3) 

The linear peak component is defined by 

Tpeak(^) z) — Hl(^) z) Tsnjooth(fc) z) , 

where the smooth component Psmooth is the linear matter power spectrum from CAMB with 
the BAO feature erased following the “sideband” method of [18]. Tabulated versions of 
Pp{k,Zrei) and T’smooth(^, Aef) are inputs in the fits. 

The key difference in our approach to BAO fitting compared to previous Lya BAO anal¬ 
yses is the treatment of the broadband distortion. Using the fe-space distortion model with 
two free parameters rather than an r-space broadband function with nine (or more) free pa¬ 
rameters, we obtain a better goodness-of-fit for the DRll auto-correlation while significantly 
reducing the number of fit parameters. 
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Figure 5. Result of applying the distortion model DC2 to the fc-weighted power spectrum Pp (left- 
hand panels) and corresponding r^-weighted correlation function (right-hand panels) at z = 2.3. 
Top panels show the results for /i = 0.9 and bottom panels for /r = 0.5. Curves indicate the undistorted 
(solid, black) and distorted (solid, blue) power spectrum and correlation function, including the effects 
of both the anisotropic non-linear broadening and the non-linear correction DNL2. The gray and blue 
bands correspond to the same parameter variations as in Fig. 4. We have adopted /3p = 1.4 and 
bp(l + (ip) = —0.351 (corresponding to bp{\ -\- (ip) = —0.336 at z = 2.25 [15]). 


A more technical difference lies in the method for evaluating the correlation function. 
The first efforts for fitting BAO in the Lya auto-correlation [16-18] and the quasar-Lya cross¬ 
correlation [21] employed a method that relied on pre-calculated templates for the I = 0,2,4 
multipoles (monopole, quadrupole and hexadecapole) of the correlation function that had 
fixed anisotropic non-linear broadening applied to the BAO peak using an approximation 
valid for (ip = 1.4. An improved method was established in [19], in which the BAO modeling 
was performed in A;-space by means of the power spectrum of the transmitted flux fraction. 
This approach allows for exact application of the anisotropic non-linear broadening. For each 
set of model parameters, the correlation function was evaluated by first expanding Pp into 
power spectrum multipoles / = 0, 2,4, and then performing a spherical Bessel transformation 
numerically of each multipole into its corresponding correlation function multipole. The non¬ 
linear broadening in principle transfers power to higher even multipoles (I = 6,8,...), but 
the contribution from these higher-order multipoles to the correlation function is negligible. 
Instead of using a multipole expansion, we employ a direct 3D FFT of Pp into the correlation 
function. This is because the distortion model transfers non-negligible power to all even 
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multipoles and thus the sum of multipole contributions to the correlation function does not 
converge for the I values {I < 14) that can be robustly evaluated in the numerical spherical 
Bessel transformation. 

The DRll auto-correlation measurement uses 50 bins of width 4 Mpc//i, spanning the 
range from 0 to 200 Mpc/h, for each coordinate ry and r±. The total number of (r||,r_L) 
bins is 2500. The fitting range is 40 < r < 180 Mpc/h, reducing the number of data bins 
in the fits to 1515. We choose 6_p(l + 13p) rather than bp to be one of the fit parameters, 
because it produces much less error correlation with /3p [15]. The six free parameters are: 
/3_p, bp{l + (3p), ay, q;_l, kc and pc- 

5 Tests on mock data sets 

In this section, we test the BAO model by fitting it to the auto-correlation functions measured 
for DRll mock data sets. Specifically, we compare the goodness-of-fit and the ability to 
recover the mock input values for the two distortion models. We start by introducing the 
mock samples that were generated for these tests. 

5.1 Mock data sets 

The method for producing mock data sets with the BOSS DRll survey geometry and sim¬ 
ulated instrumental characteristics is described in [32]. Following the method of [33], the 
absorption field for each mock sample was generated using an input Pp with I3p = 1.4 and 
bp{l + (3p) = —0.336 at z = 2.25, and a cosmology matching our fiducial cosmology, with 
the exception of assuming massless neutrinos {Q-yh? = 0). The redshift evolution of Pp was 
assumed to follow equation (2.6). Non-linear effects were added to Pp according to the model 
DNLl with the parameter values specified in section 2.2. Anisotropic non-linear broadening 
was not included for the mock samples. 

We use two types of mock samples, referred to here as standard and noiseless mocks. 
The standard mocks were produced to emulate the characteristics of the real DRll quasar 
spectra by adding instrumental noise and the effects of sky residuals, mis-calibration of the 
flux and mis-estimation of the noise from the BOSS data reduction pipeline to all the spectra 
[32]. However, absorption from intervening damped Lya (DLA) absorbers, Lyman limit 
systems (LLS) and metals were not included. A total of 100 independent realizations of the 
standard mock samples were produced. In addition, 100 noiseless mock samples using the 
same absorption field as the standard mocks but without any added instrumental noise or 
observational systematics were produced for the purpose of testing the BAO model and the 
fitting method under idealized conditions. 

Each mock sample was analyzed to measure the auto-correlation and covariance matrix 
as described in [19]. The Lya forest continua for the standard mocks were fitted using the 
“C2” method. Since the standard mocks were previously used for testing BAO fitting for 
models with the r-space broadband function in [19], we can directly compare those results 
to the results obtained here using the A:-space distortion model. For the noiseless mocks, we 
choose to use the true continua generated for each spectrum instead of performing the con¬ 
tinuum fitting. The auto-correlations for the noiseless mocks thus do not feature broadband 
distortions {Dq = 1) and are used for cross-checking the consistency with the input values 
in the absence of distortions. 
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Figure 6. Variation of the effective redshift for the standard mocks. Left: Effective redshift for each 
bin (r||,r_L). Right: Effective redshift as a function of ry for r± = 102 Mpc/h. The red solid line 
shows the best fit using equation (5.1). 


5.2 Effective redshift correction 

We adhere to the definition of the (overall) effective redshift of the BAO measurement as the 
average redshift of all the pixel pairs contributing to the correlation function in the region 
80 < r < 120 Mpc//i. The values are = 2.322 for the standard mocks and Zgg = 2.244 
for the noiseless mocks. The difference is due to the different noise and added observational 
systematics that change the weighting of the pixel pairs. Figure 6 shows the effective redshift 
for each bin in the (r||,r_L) plane for the standard mocks. The value of z^g increases in the 
ry direction, but shows negligible change in the r±_ direction. This variation along the line of 
sight is due to the increasing relative contribution from quasars at higher redshift which have 
longer forest lengths. If this effect is not corrected for in the hts, it leads to percent-level 
systematic offsets in the values oi (3p and 6i?(l + /^i?). 

We model the variation of the effective redshift using a second-order polynomial in ry, 

« (lOO Mpc/ft) + (lOO lipc/J ■ 

The right panel of figure 6 gives the effective redshift for the slice along r±_ = 102 Mpc//i. 
Equation (5.1) provides a good ht to the data points. We repeat the same exercise for the 
noiseless mocks as well as for the DRll data. The best-fit coefficients zq, z\ and Z 2 are 
summarized in table 2. 

5.3 BAO fitting of mock data sets 

We perform hts to the auto-correlation functions for 100 noiseless mock samples without 
distortion modeling {Dq = 1) and for 100 standard mock samples using distortion models 
DCl and DC2. All the hts feature the DNLl model for non-linear corrections, no anisotopic 
non-linear broadening (Sy = S_i_ = 0) and a linear matter power spectrum calculated for 
the mock cosmology at z^ei = 2.25, in accordance with the model for Pp of the mocks. 
The goal of this exercise is to recover the hducial values ay = 1, a_L = 1, f3p = 1.4 and 
bp{l + (3p) = —0.336 at z = 2.25, on average. 
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Sample 

ZQ 

Zl 

Z2 

noiseless mocks 

2.1832 

0.0221 

0.0068 

standard mocks 

2.3161 

0.0053 

0.0086 

DRll data 

2.3245 

0.0051 

0.0080 


Table 2. Best-fit coefficients for the effective redshift correction of equation (5.1) for the noiseless 
and standard mock samples and the DRll data. 


Analysis 

ay 

a± 

Pf 

6ir(l ± I3f) 

fiducial 

1.0 

1.0 

1.4 

-0.336 

noiseless mocks 

0.9993 ±0.0014 

1.0064 ±0.0023 

1.390 ±0.007 

-0.3365 ± 0.0003 

std. mocks, DCl 

1.0067 ±0.0029 

0.9984 ± 0.0054 

1.493 ±0.013 

-0.3216 ±0.0005 

std. mocks, DC2 

1.0047 ±0.0029 

1.0030 ±0.0054 

1.397 ±0.009 

-0.3347 ±0.0004 


Table 3. Summary of the results for fits to 100 noiseless mock samples and 100 standard mock 
samples. The parameter values shown are unweighted mean values with la errors and should be 
compared to the fiducial values that were used for generating the mocks. DC2 uses the fixed value 
Pc = 3, as motivated in the text. 


Starting with the results for the 100 noiseless mocks, table 3 summarizes the unweighted 
mean values with Icr errors for each parameter. While we are able to recover the correct mean 
value of ay, we notice a small but statistically significant positive offset in the mean value 
of a_L. This offset represents a systematic error that is not due to broadband modeling, 
spectral noise or the decoupling of the peak, and persists when using the fitting method 
of [19]^. There are a couple of approximations in the generation of the mocks that could 
explain the offset. First, the code to generate the mocks is based on the assumption that 
the different lines of sight are parallel. In order to add the correct geometry in our mock 
skewers, we have to generate all the skewers at different redshifts (four in the mocks used for 
this analysis), and use a linear interpolation as a function of redshift to generate a value for 
each pixel. Second, the code to generate the mocks uses the small angle approximation to 
compute the transverse separations between two pixels. The analysis code, however, uses a 
better approximation that is correct at second order in the angular separation. A sub-percent 
systematic error for a± is less than one fifth of the statistical uncertainty reported in our 
DRll analysis, and therefore can be regarded as a satisfactory performance. The recovered 
mean values of /3 f and 6 _f(1 -|- /3p) are consistent with their fiducial values. We note that 
fits with the distortion model included yield Dc —?■ 1 {kc —?• 0), as expected, but its built-in 
property of vanishing at A:|| = 0 introduces a bias in the recovered parameters. 

For the analysis of the 100 standard mocks, we begin by investigating the performances 
of the distortion models with respect to /3 f and bF{l + Pf). Figure 7 displays the best-fit 
parameter values of Pf, + Pf) and kc as a function of pc- Both distortion models 
exhibit a correlation where /3f and 6_p(l -|- Pf) increase with pc- For DC2, the degeneracy 
between kc and pc causes the fits to 52 of the mock samples to fail to converge. The 
degeneracy is evident in the bottom panel where there is a clear correlation between the 
distortion parameters. Moving forward, we fix the pc parameter for DC2 in order to remove 
this degeneracy to improve stability and performance of the fits. The value pc = 3 makes 

^This systematic error was not discovered in [19] where only standard mocks were analyzed. 
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Figure 7. Best-fit parameter values for f3F, &f(1 + I3f) and kc as a function of pc for fits to the 
standard mocks. Results are shown for 100 mock samples for DCl (red triangles) and 48 mock samples 
for DC2 (blue dots). The black dashed lines indicate the fiducial values that were used for generating 
the mocks. Errors bars have been omitted for clarity. 


both Pf and 6_p(l + Pf) agree with their fiducial values. DCl does not suffer from the same 
degeneracy, and there is no value of pc that produces agreement for both and 6^(1 + /^f)) 
so we will continue to treat it as a two-parameter model. 

The main results for the 100 standard mocks are presented in table 3. Both distortion 
models give mean values of ay and a±_ that are statistically consistent with the fiducial values 
(although also consistent with the offset value found for the noiseless mocks). It is reassuring 
to find that neither distortion model biases the BAO peak position. The corresponding values 
obtained for BAO models with the r-space broadband function in [19] are ay = 1.003 ±0.003 
and a_L = 1.008 ± 0.006, in agreement with our results. DC2 is also successful in recovering 
values of and 6^(1 ± I^f) that agree with the fiducial ones to 0.2% and 0.4%, respectively. 
We emphasize that similar results are obtained for DC2 even if pc is a free parameter (see 
figure 7). DCl, on the other hand, yields systematically higher values by about 7% for fSp 
and 4% for 5^(1 ± Pf)- We have additionally performed fits to the correlation function 
obtained from stacking all the standard mocks and find parameter results consistent with 
those in table 3. 

Figure 8 presents histograms of the minimum values for fits to the standard mocks. 
The mean values are comparable, with = 1531.6 for 1509 degrees of freedom for DCl, 
and x^ = 1534.7 for 1510 degrees of freedom for DC2. This is a sound result in support of 
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Figure 8. Distribution of the minimum values for fits to the 100 standard mocks using DCl (left) 
and DC2 with pc = 3 (right). Mean values are x^ = 1531.6 for 1509 degrees of freedom for DCl, and 
X^ = 1534.7 for 1510 degrees of freedom for DC2. 


the distortion modeling, and a significant improvement over the mean value = 1572.8 for 
1502 degrees of freedom obtained for the BAO model with the r-space broadband function 
in [19]. 

The best-fit parameter values and their la errors are displayed in figure 9. Results are 
presented for 84 mock samples for DCl and 93 mock samples for DC2, for which the fit was 
able to accurately estimate the parameter covariance matrices. The data points are scattered 
symmetrically around the hducial values for DC2, whereas DCl produces an offset for /dp 
and bp(l -|- /3p). The parameter errors are quite comparable between the models, and their 
mean values for DC2 are ajs^ = 0.081, = 0.0043, cJq|| = 0.025 and aa^ = 0.042. 

Finally, figure 10 shows scatter plots for the best-fit parameters. The data points are 
aligned horizontally for /3p and bp(l -j- /3p), indicating that these parameters are uncorre¬ 
lated. The anticorrelation between the recovered ay and a± agrees with the result of [19]. 
Comparing the values on a mock-by-mock basis, ay is typically higher for DCl by 0.1% and 
a±_ is typically higher for DC2 by 0.3%, but larger differences on the percent-level occur for 
a few mocks. 

6 Results 

Fits to the DRll auto-correlation feature the DNL2 model for non-linear corrections, hxed 
anisotropic non-linear broadening of the BAO peak and a linear matter power spectrum 
calculated for the fiducial cosmology at z^et = 2.3. The effective redshift of the measurement 
is Zes = 2.34 and its correction in the ry direction is described by equation (5.1) using the 
coefficients listed in table 2. Based on the results for the standard mocks in section 5.3, the 
baseline ht for deriving our main results uses the distortion model DC2 with pc = 3. 

Figure 11 displays the DRll auto-correlation measurement in the (rjprx) plane for the 
fitting range 40 < r < 180 Mpc//i, along with the best-fit BAO model for the baseline 
ht. Table 4 presents the best-ht results for the baseline ht along with several alternative 
hts. Parameter errors are estimated assuming the likelihood function is Gaussian at la (a 
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Figure 9. Best-fit parameter values and errors for fits to the standard mocks using DCl (red triangles) 
and DC2 with pc = 3 (blue dots). Results are shown for 84 and 93 mock samples with accurate 
parameter covariance matrices for DCl and DC2, respectively. The black dashed lines indicate the 
fiducial values that were used for generating the mocks. 



Figure 10. Plots of the best-fit parameter values for fits to the 100 standard mocks using DCl (red 
triangles) and DC2 with pQ = 3 (blue dots). The black dashed lines indicate the fiducial values that 
were used for generating the mocks. Error bars have been omitted for clarity. 


more rigorous error estimation for the baseline fit follows later in this section). The last 
row of the table (denoted “Delubaclh”) is the results for the standard fit with the r-space 
broadband function in [19] and serves as the main comparison. Both distortion models 
provide significantly better fits to the data, despite using fewer free parameters. The results 
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Figure 11. Measured DRll auto-correlation r'^^{ri\,r±) (left panel) and best-fit model using DC2 
with Pc = 3 (right panel), as a function of radial (ry) and transverse (r_L) separations. The color scale 
saturates at = 1 and = —2.5. 


Analysis 

ay 

a L 

Pf 

bpil ± Pf) 

X^IDOF 

DC2 

1.055 ±0.037 

0.963 ±0.062 

1.39 ±0.10 

-0.374 ± 0.005 

1484.2/1510 

free pc 

1.054 ±0.037 

0.963 ±0.062 

1.37 ±0.09 

-0.375 ± 0.006 

1484.2/1509 

-Dnl = 1 

1.055 ±0.036 

0.963 ±0.061 

1.33 ±0.09 

-0.377 ± 0.005 

1481.5/1510 

DCl 

1.057 ±0.032 

0.962 ± 0.059 

1.31 ±0.09 

-0.363 ± 0.007 

1484.7/1509 

Delubacl5 

1.054 ±0.031 

0.973 ±0.051 

1.50 ±0.47 

-0.402 ± 0.024 

1499.1/1502 


Table 4. Best-fit results for fits to the DRll auto-correlation. Our baseline fit uses DC2 withpc = 3 
and is presented on the first row. Delubacl5 refers to the standard fit in [19] that uses an r-space 
broadband function to accommodate the broadband distortions. 

for ay and Q!_l are consistent between all models, confirming that the measurement of the 
BAO peak position is robust to the different treatments of the broadband distortion. The 
largest benefit of using the distortion models is the much more precise constraints obtained 
for Pf and &p(l -|- /3f)- The magnitudes of the statistical error estimated for the best fit 
for these parameters are reduced by a factor of hve when using DC2. Compared to the 
best-ht values for DC2, DCl yields values that deviate by about 6% for and by 3% 
for 6p(l + Pf), compatible with the magnitude of the systematic errors estimated for this 
model in section 5.3. Two of the alternative fits vary the assumptions for the baseline ht. 
Allowing Pc to be a free parameter does not improve the goodness-of-fit or change the best-fit 
parameters notably; the reason is that the best-fit value pc = 2.57 ± 0.27 is not far from 
Pc = 3. Ignoring the non-linear correction (Dnl = 1) produces no signihcant difference 
either. 

The best-fit distortion models are presented in figure 12, with the parameter values 
specihed in the figure caption. The values of kc are larger than those found for the standard 
mocks, which suggests that there are additional sources of broadband distortion present in 
the data that are absorbed by the distortion models. 
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Figure 12. Best-fit distortion models DCl (red dashed) and DC2 with pc = 3 (blue solid) for fits 
to the DRll auto-correlation. The best-fit distortion parameters are kc = 0.030 ± 0.002 h/Mpc and 
PC = 0.847 ± 0.075 for DCl, and kc = 0.0050 ± 0.0001 h/Mpc for DC2. 
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Figure 13. Constraints on ay and aj_. Contours show 68.3% (Ax^ = 2.30), 95.45% (Ax^ = 6.18) and 
99.73% (Ax^ = 11.83) confidence contours from the DRll auto-correlation using the distortion model 
DC2 with Pc = 3. The overlayed dashed contours are for the standard fit in [19] that uses an r-space 
broadband function to accommodate the broadband distortions. Results are in good agreement. 



Figure 13 displays the 68.3% (Ay^ = 2.30), 95.45% (Ay^ = 6.18) and 99.73% (Ay^ = 
11.83 confidence contours for the scale factors ay and q;_l for the baseline fit, and compares 
them to the corresponding contours obtained in [19]. The results are in good agreement, 
although the distortion model produces slightly larger contours that are more Gaussian in 
shape. Therefore, the constraints on the Hubble parameter and the angular diameter dis¬ 
tance reported in [19], and consequently the cosmological implications, are unchanged. This 
conclusion is further illustrated in figure 14, where we combine the auto-correlation con- 
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Figure 14. Constraints on Dn/rd and DA/rd- Contours show 68.3% (Ax^ = 2.30) and 95.45% 
(Ax^ = 6.18) confidence contours from the DRll auto-correlation using the distortion model DC2 
with Pc = 3 (blue), the DRll quasar-Lya cross-correlation [21] (red), and the combined constraints 
(black). The dashed contours are the results of [19], in agreement with our constraints. The dotted 
lines indicate the fiducial values Dfj/rd = 8.708 and DA/rd = 11.59 at z = 2.34. 


tours with the contours from the quasar-Lya cross-correlation [21]. The constraints on ay 
and Q!_l have been converted into constraints on Dn/rd and DA/rd using the fiducial values 
Dn/rd = 8.708 and DA/rd = 11.59 at z = 2.34. Our combined contours are in excellent 
agreement with the result of [19]. 

Figure 15 shows the 68.3%, 95.45% and 99.73% confidence contours for I3f and bpi^+I^F) 
at z = 2.3 for the baseline fit. The marginalized 68.3% (Ax^ = 1), 95.45% (Ax^ = 4) and 
99.73% (Ax^ = 9) constraints are 


I3f = 1.39 


-1-0.11 -1-0.24 +0.38 
-0.10 -0.19 -0.28 


bF{l + Pf) = -0.374 


+0.007 +0.013 +0.020 
-0.007 -0.014 -0.022 ' 


( 6 . 1 ) 


In comparison, the marginalized 68.3% constraints for the Delubacl5 fit are /3f = 1.47^’^q™ 
and 6^(1 -b /3f) = -0.400+q;o5q, and thus the statistical errors are reduced by more than 
a factor of seven for our baseline fit. The overlayed 68.3% confidence contours in figure 15 
demonstrate how the constraints vary when changing to pc = 4 and pc = 2.5. These 
values yielded acceptable agreement for 13f and 6i?(l -|- Pf), respectively, for the standard 
mocks. This illustrates the (maximum) level of systematic error associated with the distor¬ 
tion model. Increasing pc moves Pf and 6i?(l -|- Pf) to higher values, whereas decreasing pc 
gives lower values, echoing the correlation observed in the results for the standard mocks (see 
figure 7). The value Pf = 1-4 assumed for the mock data sets and predicted in hydrodynam- 
ical simulations [27] agrees well with our measurement. Interestingly, Pf = 0.71^Q^g and 
6 ^ 7(1 + Pf) = —0.351 ± 0.013 (corresponding to 6i?(l -|- Pf) = —0.336 at z = 2.25; assumed 
for the mock data sets) reported by [15] for smaller scales 20 < r < 100 Mpc//i are in moder- 
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Figure 15. Constraints on and &f(1+/3f) at z = 2.3. Contours show 68.3% (Ax^ = 2.30), 95.45% 
(Ax^ = 6.18) and 99.73% (Ax^ = 11.83) confidence contours from the DRll auto-correlation using 
the distortion model DC2 with pc = 3. The overlayed 68.3% confidence contours are for pc = 2.5 
(dashed) and pc = 4 (dotted). Scale factors have been fixed to their best-fit values ay = 1.055 and 
a_L = 0.963. 



ate tension with our constraints. We also note that the best-fit value Pf = l-12]'iQ2g for the 
quasar-Lya cross-correlation [21] is compatible within 68% conhdence with our measurement. 

Lastly, the measured values of Pf and 6 f(1+/3f) at z = 2.3, with a correlation coefficient 
of -0.183, can be converted through error propagation into constraints on the Lya forest bias 
parameter 

=-0.157 ±0.007 , (6.2) 

and the Lya forest velocity bias parameter = bFPF/f, 

6^ =-0.226 ±0.009 , (6.3) 

where we have adopted / = 0.961 at this redshift for our fiducial cosmology. These measured 
values of and by are somewhat lower (higher in absolute value) than predictions from 
hydrodynamical simulations [13, 27] and theoretical models [34]. 

7 Discussion 

In this paper, we have described a /c-space model for the broadband distortion of the three- 
dimensional Lya forest auto-correlation function induced by quasar continuum fitting. The 
distortion model improves on the previous treatment of broadband distortions in BAO fitting 
by providing a lower x^ using fewer parameters for the fit to the BOSS DRll Lya auto¬ 
correlation. More importantly, it enables a measurement of the Lya forest redshift-space 
distortion parameter Pf to a, statistical precision of 8% and the combination with the linear 
bias parameter 6 f(1±/3f) to 2%, a reduction of the statistical errors by more than a factor 
of seven. The systematic errors on these parameters associated with the distortion model 
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are estimated using mock data sets to be less than 0.5%. Regarding the BAO result, we 
confirm that the measurement of the peak position is robust to the different treatments of 
the broadband distortion. Our results for ay and a±_ are in good agreement with [19], and 
thus the cosmological implications of the DRll Lya BAO measurement remain the same. 

Even though broadband distortions are primarily due to continuum fitting, additional 
distortions of the auto-correlation could arise from instrumental systematics (spurious corre¬ 
lations from, e.g., sky residuals and spectro-photometric calibration introduced by the data 
reduction pipeline), large-scale intensity fluctuations in the cosmic ionizing background ra¬ 
diation [35, 36], large-scale temperature fluctuations in the intergalactic medium [37], and 
absorption by high column density systems and metals [32, 38]. The fact that we obtain 
a good fit to the auto-correlation with the current BAO model suggests that these effects 
are either small, serendipitously cancel each other, or are absorbed by the distortion model. 
Regardless of the reason, it is important for future studies to better understand how these 
effects modify the auto-correlation and then correct for any significant contaminations and 
rehne the BAO model to include the relevant astrophysical effects. Improvements of the data 
reduction pipeline are anticipated for upcoming BAO analyses [39] and should alleviate or 
remove potential instrumental systematics in the measurement of the BAO peak position 
and the broadband shape. While DBAs are already accounted for in the analysis, absorber 
systems with lower column densities are difficult to reliably identify in noisy spectra and 
thus remain a potential source of systematic error. The effect of such absorbers on the 
auto-correlation is to decrease the measured values oi fip and 6^? by a few percent [38] and 
was proposed in [15] as an explanation for the discrepancy between the measured /Sp and 
theoretical predictions from numerical simulations of the Lya forest [13, 27]. The scale de¬ 
pendence introduced by fluctuations in the ionizing background leads to a suppression of the 
power spectrum on large scales. Since this result is degenerate with the continuum fitting 
distortion, it can be absorbed by the distortion model [40]. 

Given the improved modeling of the broadband shape, it is desirable to extend the fitting 
range to smaller separations (r < 40 Mpc//i) and thereby further reduce the statistical errors 
on /3f and 6 _ f (1 + Pf)- Such a fit finds lower [3f and higher 5i?(l -|- /3i?), changes in the 
directions of the values reported in [15]. However, the distortion models considered in our 
analysis do not give unbiased measurements for the larger fitting range. The systematic 
errors on /3f and 6 _f(1 + /3f) increase to 2-3% for fits with r > 20 Mpc//i, roughly the same 
magnitude as the statistical errors, and thus a more careful study involving modifications 
to the distortion model parameterizations is needed to improve on this result. Furthermore, 
the use of smaller scales would require more attention to the non-linear modeling. The 3D 
EFT resolution employed in this analysis restricts the wavenumber range to A: < 0.79 h/Mpc, 
which omits the contribution from the small-scale power spectrum in the calculation of the 
correlation function. Another interesting possibility that requires further refining of the 
broadband modeling is to constrain the scale factors ay and a±_ using the full shape of 
the correlation function in the Alcock-Paczyhski test [41]. Fits without decoupling of the 
peak from the broadband component could potentially reduce the statistical errors on the 
measured scale factors by a factor of two. 

We expect the general approach of the distortion model to also apply for the quasar- 
Lya cross-correlation, albeit with some modifications. Since the length of the Lya forest is a 
function of redshift, and the contributions to the cross-correlation function for positive radial 
separations (absorption behind the quasar) are produced by spectra with predominantly 
higher redshifts compared to negative radial separations, the broadband distortion becomes 
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asymmetric with respect to ry = 0 (see figure 2 of [21]). We therefore anticipate that the 
distortion model needs to be generalized to accommodate this asymmetry, in order to achieve 
optimal resnlts for BAO fitting of the cross-correlation. 

The BOSS data set has opened the possibility to measnre the expansion rate and ge¬ 
ometry of the Universe at redshift z > 2 using BAO in the Lya forest. The precision of 
the BAO measnrement, the level of sophistication in the fitting method and our understand¬ 
ing of the systematic errors have improved for each analysis of the growing data set, and 
the final BOSS data release [12] (DR12) will provide the best measurement yet of the BAO 
scale at z ~ 2.3. The advances in broadband distortion modeling presented in this paper 
now also enable significantly more accurate measurements of the large-scale linear bias and 
redshift-space distortion of the Lya forest. 
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